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^ ABSTRACT 



Remote observing of exoplanetary atmospheres is now possible, offering us 
access to circulation regimes unlike any of the familiar Solar System cases. At- 
mospheric circulation models are being developed to study these new regimes 
but model validations and intercomparisons are needed to establish their consis- 
tency and accuracy. To this end, we present a simple Earth-like validation of 
the pseudo-spectral solver of meteorological equations called IGCM (Intermedi- 
ate General Circulation Model), based on Newtonian relaxation to a prescribed 
^ latitudinal profile of equilibrium temperatures. We then describe a straightfor- 

ward and idealized model extension to the atmospheric flow on a hot Jupiter 
<N with the same IGCM solver. This shallow, three-dimensional hot Jupiter model 

I is based on Newtonian relaxation to a permanent day-night pattern of equilib- 



rium temperatures and the absence of surface drag. The baroclinic regim^ of 
the Earth's lower atmosphere is contrasted with the more barotropic regime of 
0\ the simulated hot Jupiter flow. For plausible conditions at the 0.1-1 bar pres- 

QQ sure level on HD 209458b, the simulated flow is characterized by unsteadiness, 

O subsonic wind speeds, a zonally-perturbed superrotating equatorial jet and large 

scale polar vortices. Violation of the Rayleigh-Kuo inflexion point criterion on 
^ the flanks of the accelerating equatorial jet indicates that barotropic (horizontal 

shear) instabilities may be important dynamical features of the simulated flow. 
Similarities and differences with previously published simulated hot Jupiter flows 
are briefly noted. 

Subject headings: Stars: Planetary Systems, Stars: Atmospheres, Turbulence, 
Infrared: General, Infrared: Stars 



"'^Contours of pressure and density are misaligned in a baroclinic flow, while they are aligned in a barotropic 
flow. Strongly baroclinic flows are susceptible to baroclinic instabilities which, in essence, are slanted versions 
of convection (see the review by Showman et al.|[2008 ) . 
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Introduction 



Since tlie first discoveries of extrasolar giant planets around nearby sun-like stars with 
Doppler velocimetry (Mayor Sz Queloz 1995; Marcy & Butler 1996), extrasolar planet re- 
search has experienced a spectacular series of observational breakthroughs. In recent years, 
progress has been particularly rapid with a subset of close-in extrasolar planets found to 



transit the disk of their host star (e.g., Charbonneau 2008 [Deming 2008). 



The first transit observation of a hot Jupiter (Charbonneau et al. 2000; Henry et al. 
2000) has been followed by many such measurements, including a growing number of new 



close-in planet discoveries based on transit searches (e.g., Torres et al. 2008 ). For the brightest 
nearby systems with transiting planets, this has also enabled the detection of the planetary 
thermal flux occulted at secondary eclipse (e.g., Deming et al. 2005, 2006; Charbonneau 
et al. 2005; Harrington et al. 2007), the detection of day/night temperature variations 
through IR phase curve monitoring (e.g. Harrington et al. 2006; Knutson et al. 2007, 
2008b; Cowan et al. 2007), IR spectral measurements (Grillmair et al. 2007; Richardson 
et al. 2007; Knutson et al. 2008a) as well as a variety of additional constraints based on 
transit spectroscopic studies (Tinetti et al. 2007; Ehrenreich et al. 2007; Swain et al. 2008; 
Barman 2007; Redfield et al. 2008; Pont et al. 2008). Overall, these results have contributed 
to a shift in focus from the detection of extrasolar planets to the characterization of their 



physical atmospheric properties (e.g., ' 


meager et al. 


2005 


Burrows et al. 


2006 


Marley et al. 


2007 


Fortney & Marley 


2007 


Barman 


2008 


Burrows et al. 


2008 


)■ 



The need to interpret these astronomical data reliably, so as to infer the physical con- 
ditions present in distant exoplanetary atmospheres, has also fueled a growing atmospheric 
modeling effort. While plane-parallel radiative models have been the tools of choice to in- 
terpret these data until now, they are likely insufficient in the case of hot Jupiter/Neptune 
atmospheres. Indeed, due to their short orbital separations, all these exoplanets are expected 
be tidally-locked to their parent star (or, in some cases, pseudo-synchronized) . Tidally- locked 
planets are subject to an unusual situation of permanent, asymmetric day/night radiative 
forcing, leading to heat redistribution by atmospheric motions. Current data on hot Jupiters 
already indicate that existing one-dimensional radiative transfer models fail to capture the 
multi-dimensional nature of this atmospheric regime (e.g., Seager et al. 2005; Knutson et 
al. 2007; 2008b; Fortney et al. 2006) and, when included, the effects of atmospheric heat 



redistribution are usually accounted for only indirectly (e.g., Iro et al. 2005; Seager et al. 



2005 Barman 2008 Burrows et al. 2008 Fortney et al. 2008). To adequately address this 



circulation regime and interpret the growing data set, multi-dimensional, coupled radiation- 



hydrodynamics models of these atmospheres are required (see the review by Showman et al. 
2^081). 
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Hot Jupiters have high atmospheric temperatures, slow rotation rates and an unusually 
permanent pattern of asymmetric day/night radiative forcing. Recognizing that this consti- 
tutes an interesting problem in atmospheric dynamics, several groups have explored plausible 
circulation regimes on these planets using different modeling approaches and assumptions 



(Showman & Guillot 2002 


Cho et al. 


2003 


Vfenou et al. 2003; Cooper & Showman 2005 


2006 Langton & Laughlin 


2007 


Cho et al. 


2008 Dobbs-Dixon & Lin| 2008 Langton & 


Laughlin 2008; Showman et al. 


2008) 


. These investigations have exhibited significant di- 



versify in flow results, for which there currently is no simple unifying explanation (see the 



discussion in Showman et al. 2008 ) 



Atmospheres in motion are non-linear, presumably turbulent flows. In addition, radia- 
tively active species in an atmosphere are advected by the flow and in doing so non-linearly 
couple the flow to its effective source of radiative forcing. Simulated atmospheric flows are 
thus generally expected to be sensitive to various numerical and physical details of any spe- 
cific implementation. In the exoplanetary context, with rather limited direct information on 
the atmospheres studied, this means that simple parameterized models isolating key features 
of the simulated flow are important in helping us understand the general behavior of remote 
atmospheric flows. With this in mind, we present here idealized atmospheric models which 
emphasize dynamical aspects under linear forcing conditions. These simplified models do 
not address any detailed aspect of the radiative or chemical structure of the atmospheres 
under consideration. 

In general, different numerical implementations of a specific atmospheric problem will 
not lead to identical results. For this reason, model validations and inter-comparisons are 
standard practices in atmospheric science. This is particularly true of the complex hydro- 



dynamic solvers known as dynamical cores (e.g.. Held & Suarez 1994). The development 
of dynamical cores to solve the equations satisfied by an hydrostatic atmosphere in motion 
has been a major enterprise. Dynamical cores currently used in meteorological and climate 
models are the results of years of refinements to guarantee stability, efficiency, and accurate 
conservation of mass, momentum and energy (e.g., Hoskins & Simmons 1975). In the exo- 



planetary context, different modeling groups have used various hydrodynamic solvers. Some 
are new and httle tested while others are old, well-tested but only validated under conditions 
appropriate for Solar System planetary atmospheres. As a result, intercomparisons and val- 
idations of dynamical cores on idealized atmospheric problems tailored for exoplanets will 
be important in the future to assess, both qualitatively and quantitatively, how reliable our 
interpretations of exoplanetary data can be. 

As a first step in this direction, we present here a simple Earth-like validation of the 
IGCM dynamical core (in §3), followed by a direct extension of this model to the atmospheric 
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flow on a hot Jupiter (in §4). 



2. The IGCM Dynamical Core 



The IGCM dynamical core has its origins in the development of accurate numerical 



solvers for meteorological studies (Hoskins & Simmons 1975). It is a well tested and docu- 



mented solver, which has been used for studies as diverse as the climate of the Earth (Forster 



et al. 2000), the meteorology of Mars (Joshi et al. 1995; Collins & James 1995) and the at- 



mospheric circulation regime of tidally-locked terrestrial exoplanets (Joshi et al. 1997). We 
summarize the main features of this dynamical solver here. 



2.1. Equations of Motion 



We consider a coordinate system rotating with the same angular velocity, Qp, as the 
planet. Longitude A ranges from to 2tt in the eastward direction and latitude ranges from 
—tt/2 at the south pole to +7r/2 at the north pole. From the general fluid equations satisfied 
by an inviscid ideal gas, one derives the primitive equations of meteorology in the traditional 
approximation, which are satisfied on large horizontal scale^ by a hydrostatic atmosphere 
in motion. Using pressure, p, as the vertical coordinate, equations for the conservation of 
momentum, mass and energy can be written, respectively. 



dv u tan (f) 



dt 



kx V 

dp 
duj 
dp 
dT 



-Vp<l>-/kxv, 

1 
P 



-Vp ■ V , 



UJ 



+ Q, 



(la) 
(lb) 
(Ic) 
(Id) 



where the Lagrangian derivative is 

d 
dt 



— V V 

dt ^ dp' 



(le) 



^The primitive equations cease to be satisfied by motions with horizontal scales approaching the atmo- 
spheric vertical scale height (e.g., Holton 1992; Vallis 2006). 
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V = [u, v) is the (eastward, northward) horizontal velocity in the rotating frame, Rp is the 
planetary radius and u = dp/dt is the vertical "velocity." The geopotential $ = gz, where z 
is the height above a fiducial pressure surface and g is the gravitational acceleration, which is 
assumed to be constant and includes the contribution from the centrifugal acceleration. The 
unit vector k is normal to the "planetary surface." The Coriolis parameter / = 2f2psin0 
is the projection of the planetary vorticity vector, 2f2p, onto k. The operator Vp is the 
horizontal gradient on a constant-j9 surface. In the energy equation, Cp is the heat capacity 
of the atmospheric gas at constant pressure and Q is the diabatic heating/cooling rate. The 
above equations are closed with the ideal gas law, p = pTZT, as the equation of state, where 
TZ is the atmospheric perfect gas constant. The validity of the primitive equations in the 
traditional approximation for the study of hot Jupiter atmospheric fiows has been discussed. 



for instance, by Showman et al. (2008) and Cho et al. (2008). 



The IGCM dynamical core solves these equations in a form that involves the vorticity 
and divergence of the horizontal fiow. Using scaled units of length (oc Rp), tim^ (oc fip^) 
and temperature (cx R^Vt^/lZ), as well as rescaled zonal and meridional velocities 

U = MCOS0, 

V = V cos 0, 

the horizontal momentum equations specifically solved are the equation for the vertical 
component of absolute vorticity, (, 

dC _ 1 dJ^v _ dTu 

dt 1 — /i^ d\ dfi 
and the equation for the horizontal divergence, D, 



~t~ Qvor ~t~ Q 



hyp 
vor ' 



(2) 



dD 
~dt 



dJ^u dJ^v 

+ 



1 — jj? dX djj, 



2(1 



$ + Tref Inpsurf 



+ Qdiv + Q5L^ (3) 



where 



c 

D 



sm (p 
2/i + 



dv 



1 



ducoscf) 



Rp cos d\ Rp cos d(j) 



du 



+ 



dvcosc 



Rp cos dX Rp cos 
dU ^ 9 In »surf 



-f/C - 0-- Ta(l - /i ^ 



dV 



^All times are subsequently quoted in units of a planetary day equal to 2Tr /flp. 
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Conservation of energy is expressed via the temperature equation 

^ = -T^2&,UT^ - ^VT^ + DT^ -a^ + n^ + Qr + , (4) 
ot 1 — /i^ OA ofx da p 

where the temperature T = T^ef + Ta is decomposed into an arbitrary (constant) reference 
value, Tref(cr), around which the anomaly, Ta(cr), is calculated. Conservation of mass is 
expressed via an equation for the vertically integrated "surface" pressure, 

dt ~ 1-/^2 d\ dfi da ^ ^ 

where psurf is the pressure at the bottom model level, due to the entire weight of the overlaying 
atmosphere. Momentum balance in the vertical satisfies the same hydrostatic equation as 
before, 

7^ = (6) 
dm a 

In all the above equations, the vertical pressure coordinate a = p/psmi is used (0 < a < 1). 
These five equations are the exact same ones as originally solved by Hoskins fc Simmons] 
(1975), except for the various additional Q terms entering the vorticity, divergence and 
temperature equations. 



2.2. Forcing, Friction and Hyperdissipation 

Our IGCM implementation contains linear prescriptions for temperature relaxation as 
well as drag on the vorticity and divergence fields. Newtonian relaxation of the temperature 
field to a prescribed equilibrium temperature profile, Teq, on a characteristic radiative time 
Trad, corresponds to 

T -T 

Qt = -^^^ (7) 

''"rad 

in equation (Q. This linear temperature relaxation scheme is a very simplified form of 
atmospheric forcing which encapsulates all the unspecified chemical-radiative physics in the 
modeled atmosphere. Rayleigh drag on the horizontal fiow field, on a characteristic friction 
time Tfric, corresponds to 



C-2/i 



Qdiv = -— (9) 

^fric 

in equations ^ and (|3]), respectively. (Note that friction is only applied to the relative 
vorticity field.) 
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In this work, Rayleigh drag is used as a simple parameterization of surface friction in 
the Earth-hke model only, while Newtonian relaxation is used as the unique form of forcing 
in all models. The vertical profile of equilibrium temperature used has two regions, a lower 
atmosphere (troposphere) where temperature decreases with height at a fixed rate, and 
an upper atmosphere (stratosphere) which is vertically isothermal. The specific relaxation 
temperature profile adopted in all the models has the form 



T:^^\z) = Ts^rf - r,,op(^stra + ^^) + J (^h,,,,[z - Z,,,,]^ + 5T,i„ (10) 



where z is the height above the lowest "surface" level, Tgnrf is the "surface" temperature, 
Ftrop is the lapse rate (= —dT/dz ) in the tropospheric region, ^stra is the height at which 
one enters the stratospheric region and STstra. is a temperature offset used to smooth out the 
transition between the finite tropospheric lapse rate and the isothermal stratosphere. 

Differential heating and cooling of atmospheric regions around the planet drives atmo- 



spheric motions. In addition to the vertical dependence in equation (10), one must also 
specify the latitudinal and longitudinal dependence of the temperature relaxation profile. In 
our models, the three-dimensional equilibrium temperature is given by 

Teq(a, A, 0) = T;-H^) + Arop(a) ATe(A, 4>), (11) 

where 

Arop(fT) = sm — (12) 

2(1 - (Tstraj 

is a height-dependent damping factor applied in the troposphere (cr > cr[-2stra]) to gradually 
reduce the latitudinal and longitudinal temperature differential over the vertical extent of 
the troposphere. In the stratosphere, this horizontal temperature differential is set to zero: 
Arop(o" < o"[^stra]) = 0. The Correspondence between z and a levels is obtained by vertical 
integration of the hydrostatic balance equation (Eq. [6]). The latitudinal and longitudinal 
dependence of the relaxation temperature profiles, ATq^X, 0), are specified separately for the 
Earth-like and hot Jupiter models in § [3] and § |4] below. 

When using the primitive equations, one assumes that effects on small scales, which 
are not adequately described by the equations, nor resolved numerically, can be parameter- 
ized within the framework of the large-scale dynamics. In particular, in a turbulent two- 
dimensional flow, enstrophjj^ is known to cascade to small scales, where three-dimensional 

'^Enstrophy is the "vortical energy" of the flow, defined by jCroi' where Croi — Crci k — Vxv is the relative 
vorticity (parallel to the vertical direction k). 
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effects eventually become important (e.g., Pedlosky 1987). In the absence of explicit dissi- 



pation in a numerical atmospheric model, enstrophy accumulates at the smallest resolved 
scales, in a process known as spectral blocking, with ever growing numerical errors. It is 
thus standard practice in atmospheric science to introduce hyperdissipation terms in the 



equations solved to alleviate the spectral blocking problem (e.g., Stephenson 1994). 



In the IGCM, hyperdissipation terms are introduced in the vorticity, divergence and 



temperature equations. 



Qhyp 



-^^diss(-l)''""^V2^°-^ (C - 2^i) 

72Nr 



-^^diss(-l)''°^^V^^^°^^Ta, 



72Nr 



(13) 
(14) 
(15) 



where A^del is an integer. (Note that hyperdissipation is only applied to the relative vorticity, 
Crei = C ~ 2/i.) Hyper-laplacian operators guarantee that only scales close to the smallest 
resolved features are selectively chosen for diffusion. The case A^del = 1, which would 
correspond to a regular diffusion operator, is considered to be too widely dissipative for 
what are effectively inviscid planetary atmospheres. 

It should be noted that hyperdissipation is a numerical tool with intrinsic limitations. 
It is not supported by any fundamental theory and it may well miss some of the "reverse" 
interactions occurring between small and large horizontal scales in a turbulent atmospheric 
flow. There is no rigorous way to chose the magnitude or the order of hyperdissipation in a 



given model (e.g., McVean 1983) and different choices have typically been made by different 



modehng groups in Earth atmospheric studies (e.g., Stephenson 1994). In the context of 
exoplanet atmospheric modeling, with a priori unknown circulation regimes, hyperdissipation 
choices should probably be considered as important free parameters of the models. In keeping 



with standard practice (McVean 1983 Stephenson 1994), we use a highly scale-selective 
-^DEL = 4 scheme and a dissipation rate adjusted with horizontal resolution so that structures 
on the smallest resolved scales are dissipated in a fraction (~ 0.25 - 0.01) of a planetary day. 



2.3. Numerical Solutions 



The IGCM solves equations with the semi-implicit pseudo-spectral method 

( Hoskins fc Simmons|1975 ). The equations are solved in spectral space horizontally and with 
a finite difference scheme in the vertical. A standard decomposition in spherical harmonics 



is used with triangular mode truncation (Orszag 1970 Eliassen et al. 1970 Washington & 



Parkinson 1995). The pseudo-spectral method is well adapted to problems with spherical 



geometry and thus does not require any special treatment at the poles. The vertical finite 
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difference scheme is based on the formulation of Simmons & Burridge (1981), which con- 



serves angular momentum and energy exactly. As is usually the case with this formulation, 
total atmospheric mass is not conserved by the vertical scheme but is instead corrected for 
at each time-step. Vertical levels are linearly spaced in a. In what follows, model resolutions 
are reported as TpLg, which corresponds to 3p + 2 spectral modes in longitude, half that in 
latitude, and q vertical levels. Time integration is performed with a semi-implicit leapfrog 



scheme, followed by a Robert-Asselin filter to control time splitting (Hoskins &; Simmons 



1975 Forster et al. 2000). A significant advantage of the semi-implicit method is that it 



solves gravity wave propagation in an implicit manner. Larger time steps are thus possible 
since CFL stability is determined by advection, rather than the speed of gravity waves. 

Boundary conditions are applied only at the top [a = 0) and bottom (a = 1) of the 
numerical domain, where da/dt = is imposed. By default, this corresponds to free-slip 
boundary conditions for the horizontal flow. All our models are started at rest, with a 
temperature profile satisfying everywhere the equilibrium profile 
A small amount of noise is introduced to break flow symmetries. 



vert / 
eq y 



in equation (10). 



2.4. Flow Representations 

There are various ways to represent a three-dimensional atmospheric flow. In the inter- 
est of clarity and conciseness, we only show cylindrical maps and zonally- averaged contour 
plots. Our temperature and velocity maps use Miller cylindrical projections, centered on the 
equator (and the substellar point, at [A, 0] = [0,0], in the case of the hot Jupiter model). 
In some cases, only a subset of all velocity vectors are shown, for better rendering. We 
also present zonally-averaged contours of zonal velocity, u, and temperature, T. These were 
obtained by performing the longitudinal average 

[X{t)] = [X{a, 0, t)] = ^l ^(^' 0' dX, (16) 

where X = u or T. 



3. Earth-Like Model 



Our Earth-like model comes as a default implementation of the IGCM solver. It has 
many similarities with the control run described by James & Gray (1986). To capture in a 
simple way Earth's annual mean conditions, a meridional gradient of relaxation temperatures 
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is imposed, such that 

ATe{X, 0) = ATe(0) = - sin^ 0) x ATep, (17) 
where the equator-pole temperature difference is set to AT^p = 60 K. Parameters for the 



vertical relaxation profile, T^^ (z) in equation (|10|), are chosen to capture approximately 



Earth's typical radiative- convective conditions. For simplicity, a single value for the radiative 
relaxation time. Trad = 15 planet days, is adopted throughout the atmosphere. Rayleigh drag 
with a friction time rfnc = 1 planet day is applied to the lowest model layer only. There 
is no account of seasonal variations or topography in the model. Stratospheric conditions 
are rather poorly captured. The complete list of parameters of our Earth-like model with 
T42L15 resolution is provided in Table [T] 



This model also share strong similarities with the classic benchmark calculation of Held 



& Suarez (1994) for dynamical core validations. A significant difference with the Held-Suarez 
benchmark is that the variations of Trad and Tfnc values with location in the atmosphere in 
that model are reduced to single values in the simpler model presented here. Our Earth-like 
model may thus be considered as a crude version of the Held-Suarez benchmark, with a 
reduced number of free parameters. 

The main features of Earth's general atmospheric circulation and the ability of idealized 
models such as the Held-Suarez benchmark to reproduce them with reasonable accuracy are 



well known (e.g.. Held & Suarez 1994). The key features of the tropospheric circulation, at 
p ^ 0.1 bar (cr ^ 0.1) levels, are: (i) an axisymmetric meridional circulation via Hadley cells 
extended from the equator to approximately ±30° in latitude, which efficiently reduce merid- 
ional temperature gradients in that region, (ii) strong baroclinic activity at mid-latitudes, 
with typically ~ 6 large scale barochnic eddies spread in longitude along a "storm track," 
(iii) a zonal circulation at the surface which is characterized by easterly (westward) trade 
winds in subtropical latitudes, westerlies (eastward wind) in mid-latitudes and weak easter- 
lies near the poles, (iv) a westerly component of the zonal wind increasing with height at all 
latitudes, until localized wind maxima known as the jet streams are reached at the ~ 0.2 bar 
(a ~ 0.2) level, with peak winds at ±40-50° latitude. 

Figures [T] and [2] illustrate how our simple Earth-like model qualitatively reproduces this 
general circulation regime. While time averages over hundreds of planet days are tradition- 



ally used to characterize the flow in the Held-Suarez benchmark (Held & Suarez 1994), we 
have chosen to present snapshots of the atmospheric flow for simplicity and to facilitate 
comparisons with the hot Jupiter model presented below. Figure [T] shows temperature and 
velocity maps at planet day 150 in the Earth-like model. The top panel exemplifies the strong 
baroclinic activity that characterizes mid-latitudes, shown here in the bottom model layer 
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at the o = 0.97 0.97 bar) level. The bottom panel exemplifies the reduced baroclinic 
activity and the formation of jet streams that characterize the upper troposphere, shown 
here at the a = 0.37 (~ 0.37 bar) level. 



Figure |2] shows zonally-averaged contours of zonal wind speeds {[u] in m s~^; top panel) 
and temperature ([T] in K; bottom panel) for the same Earth-like fiow at planet day 150 
as shown in Figure [l| While the [u] contours reproduce the main qualitative features of 



the Held-Suarez benchmark calculation (see Fig. 2 in Held & Suarez 1994), quantitative 
discrepancies emerge, most notably in the extremas of wind speeds and in the detailed 
shape of the zonal wind structure (with jet stream cores incorrectly pushed against the top 
layer in our model). We attribute these differences to the simpler nature of Newtonian 
forcing and Rayleigh drag in our model, different relaxation profiles and our focus on a 
fiow snapshot rather long-term averages. Similarly, a comparison of the [T] contours shown 



in the lower panel of Figure [2j with the corresponding Figure Ic in Held & Suarez ( 1994 ) 
reveals broad quahtative agreement (e.g., fiattened equatorial contours) but also quantitative 
discrepancies. These discrepancies can be partly attributed to the different vertical profiles 



of relaxation temperature adopted here and in Held & Suarez (1994). 



Rather than focusing on a strict reproduction of the Held-Suarez benchmark results, 
which is of limited interest for a well-tested code like the IGCM solver (see, e.g., the clima- 



tology of Forster et al. 2000), our simple Earth-like model may offer interesting insight into 
important issues of model parameterizations and target accuracies for exoplanet atmospheric 
circulation studies. By comparison with a somewhat higher complexity model like the Held- 
Suarez benchmark calculation, it provides a measure of the ability of strongly parameterized 
models to successfully reproduce the main qualitative features of an atmospheric circulation 
regime like that on Earth. We note, however, that in both our model and the Held-Suarez 
benchmark, parameters were adjusted a posteriori to reproduce a known circulation regime. 
By contrast, circulation regimes are a priori unknown on exoplanets and it may be difficult 
to determine from first principles the temperature relaxation profiles and relaxation times 
needed to adequately drive or drag the flow on a remote planet. On the other hand, it could 
also be that qualitative agreement at a level comparable to that achieved by our simple 
Earth-like model turns out to be sufficient to interpret with confidence typical remote as- 
tronomical observations of exoplanets. The issue of target accuracies for reliable exoplanet 
data interpretation has received little attention until now. We will simply note here that, 
in addition to validations and inter-comparisons, parameter-space explorations with sim- 
ple atmospheric circulation models may be important ingredients of an effective strategy to 
address this data interpretation challenge. 
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4. Shallow Hot Jupiter Model 



To capture the permanent day-night forcing conditions present on a tidally-locked hot 
Jupiter, our IGCM solver has been modified to permit horizontal gradients of relaxation 
temperatures of the form 



COS A cos </) X ATep, 



which places the substellar point at (A,0) = (0,0). In the specific hot Jupiter model pre- 
sented here, the equator-pole temperature difference is set to ATep = 300 K, which cor- 
responds to a full day-night temperature difference of 2ATep = 600 K. The amplitude of 
this day-night temperature forcing, which is an important free-parameter of our model, is 
comparable to the corresponding forcing amplitude at the ~ 1 bar level in the circulation 



model of Cooper & Showman (2005). A relatively steep, linear dependence of the equilib- 
rium relaxation temperature with the cosine of the angle away from the substellar point has 
been adopted because it accounts for the extra atmospheric depth crossed by radiation at 
inclined angles and is broadly consistent with detailed radiative transfer calculations (e.g.. 



Showman et al. 2008 ) . 



Parameters for the vertical relaxation profile, 



vert / 
eq \ 



in equation (10), are chosen to 
match approximately the profile of Iro et al. (2005) in the ~ 0.1-1 bar region for HD 209458b, 
with a constant lapse rate Ftrop = 2 x 10~"^ K m~^ and no stratosphere. A significant feature 
of our vertical relaxation profile is that the day-night temperature differential asymptots to 
zero in the uppermost modeled layers, like it does in the Earth-like model (see Eqs. 



11 



12]). 

While this choice facilitates a direct comparison of circulation regimes between the Earth- 
like model (with meridional forcing) and the hot Jupiter model (with hemispheric forcing), 
it may also be a poor assumption for a hot Jupiter atmosphere. More realistically, the day- 
night temperature differential would extend to layers higher up in the atmosphere (beyond 
those modeled), with the possible existence of a stratosphere depending on the presence of 



an absorbing compound such as TiO/VO or Sulfur (] 


Jubeny et al.||2003 


Fortney et al. 


2008 


Spiegel et al. 2009| Zahnle et al. 


2009 


)• 



We adopt a single value for the radiative relaxation time. Trad 
1.5 X 10^ 



0.5 planet day ~ 
([20051 . While the 



s, to match the radiative timescale at 1 bar from Iro et al 
radiative times are expected to vary substantially with depth in hot Jupiter atmospheres, a 
single Trad value, like in our Earth-like model, is the simplest acceptable form of radiative 
forcing in a shallow atmospheric model such as ours. No Rayleigh drag is implemented. All 
the other parameters are chosen appropriately for the hot Jupiter HD 209458b. The complete 
list of parameters of this idealized hot Jupiter model with T42L15 resolution is provided in 
Table [TJ We emphasize that the model is quite shallow in the sense that only 1-2 vertical 
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levels are present above the a = 0.1 level (given the linear-a grid and vertical resolutions 
used) and that absolutely no account is made of deeper atmospheric layers present below the 
1 bar pressure level (which bounds our model at o" = 1). Despite its great simplicity, a clear 
advantage of this shallow hot Jupiter model is that it is a direct extension of our Earth-like 
model setup to the case of a hot Jupiter atmospheric flow and thus permits straightforward 
comparisons between the two simulated circulation regimes. 

Figure [3] shows temperature and velocity maps at planet day 100 in the hot Jupiter 
model, at the o = 0.7 (top panel) and 0.37 (bottom panel) levels. The temperature fields 
and particularly the velocity fields share strong similarities at these two levels. This vertical 
flow alignment, which persists throughout the various modeled layers, together with the 
lack of any identifiable baroclinic eddies, is characteristic of a barotropic flow regime (e.g., 
Cho et al. 2003 Menou et al.||2003 Cho et al. 2008). The flow is characterized by a zonally- 



perturbed superrotating equatorial wind, flanked by dynamically active vortices, counter-jets 
at mid-latitudes and the presence of large scale polar vortices. Advection of heat away from 
the substellar point, at (A, 0) = (0,0), occurs both eastward in the equatorial regions and 
westward in mid-latitudes, where the counter-jets are present. 

Figure |4] offers another view of the circulation regime in this shallow hot Jupiter model 
with zonal averages of the zonal wind velocity ([«]) at planet day 100. Over nearly the entire 
vertical extent of the model layer, the flow exhibits the broad super-rotating (eastward) 
equatorial wind and slower westward counter-jets at mid-latitudes. Maximum wind speeds 
are ^ 1800 m s~^ and zonal averages are ^ 1300 m s~^. These values are well below 
corresponding sound speeds, which range from 2.4 to 3.1 km s~^ from top to bottom of the 
modeled region. 

Various diagnostics can be used to verify that the flow has reached a stationary state 
with respect to the imposed forcing. We have performed extended runs for up to several 
hundred planetary days and have found stationary conditions for the shallow hot Jupiter 
model under consideration. To illustrate this. Figure |5] shows the time evolution over 100 
planet days of representative velocities and temperatures at various horizontal locations on 
the a = 0.5 model level. The zonal average and maximum values of the zonal velocity u along 
the equator are shown in the top panel as solid and dashed lines, respectively. After a rapid 
acceleration phase lasting ~ 5-10 planet days, a flow stationary state is reached at planet 
day ~ 20, with signiflcant fluctuations. Throughout this evolution, flow velocities remain 
subsonic. In the bottom panel of Figure [5} the corresponding evolution of temperatures 
is shown at the sub- and antistellar points (top and bottom solid lines, respectively), at 
the east and west equatorial limbs (top and bottom dotted hues, respectively) and at the 
north and south poles (two dashed lines). Despite eastward heat advection at the equator. 
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which results in comparable temperatures at the east equatorial limb and the substellar 
point, temperatures are far from being horizontally homogeneous. Note in particular that 
temperatures around the planetary limb (dashed and dotted lines) represent a diverse range 
of physical conditions even on a fixed (cr = 0.5) pressure level (see also Fig. Isl). 



We have found that these results are broadly confirmed at different, and in particular 
higher, model resolutions. Figure |6] presents a specific test of numerical convergence for our 
results. Temperature and velocity maps at planet day 100 are shown for the same shallow 
hot Jupiter model as before, except that reduced and enhanced numerical resolutions were 
used, both horizontally and vertically. The top panel shows a map at the a = 0.5 level in 
a T21L5 model. The bottom panel shows a corresponding map, at the a = 0.52 level, in a 
T170L20 model. Values of the hyperdissipation coefficient Udiss were adjusted to 7.3 x 10^^ 
and 4.7 X 10^^ m^ s~^ in these T21L5 and T170L20 models, respectively. The overall similarity 
of these temperature and velocity maps, as well as other flow attributes (e.g., zonal wind 
contours), indicates that good convergence is already achieved around T31L10 to T42L15 
resolution. Even the T21L5 flow shares many of the global attributes of higher resolution 
simulated flows. 

As mentioned earlier, the zonally-perturbed equatorial wind and its flank vortices are 
dynamical features of the simulated flow. Figure [7] illustrates the flow unsteadiness with 
two successive temperature and velocity maps at planet days 97 and 98 in a T85L20 version 
of our shallow hot Jupiter model (with z/diss adjusted to 5.9 x 10^^ m^ s~^). The maps 
are shown at the a = 0.52 level. While this particular example was chosen to exhibit 
clear temperature and flow field variability over one planet day, dynamical variability is a 
general property of the simulated flow (see fluctuations in Fig. [s]). Nevertheless, we find 
that the displaced anticyclonic polar vortices that emerge in this hot Jupiter model do not 
experience systematic longitudinal translations, nor are they close to geostrophic balance. 



like the cyclonic circumpolar vortices discussed by Cho et al. (2003, 2008). Instead, the polar 
vortices in the present model appear to be strongly tied to the imposed day-night forcing 
and show only limited excursions away from their preferred night-side location, somewhat 
eastward of the anti-stellar point (see Figs. [3| |6]and[7]). 

The wind acceleration episode apparent in Figure [5] must be important in determining 
the nature of the stationary regime eventually achieved by the forced flow. We have found 
evidence that barotropic instabilities play a role in shaping the dynamical nature of this flow. 
During the first few planet days in our hot Jupiter model, we observe the acceleration of a 
zonal eastward wind in the equatorial regions and westward counter-jets at mid-latitudes. 
This situation is very reminiscent of the counter-acceleration caused by meridional Rossby 



wave transport that occurs in the idealized momentum-forced flow discussed by Cho et al. 
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(2008) in the context of the equivalent-barotropic formulation (see their Fig. 17). Figure |8] 
shows a temperature and velocity map at planet day 5 and level a = 0.9 in our T42L15 hot 
Jupiter model. In contrast with previous maps, this one zooms in a specific region along 
the equator, restricted +20 to +160° in longitude and ±40° in latitude. In this region, the 
flow exhibits strong horizontal (zonal) shear, as well as small scale velocity and temperature 
disturbances along the leading edge of the equatorial jet. Subsequently, we observe a rapid 
thinning of the jet, a breaking of equatorial symmetry by planet day 6-7 and the emergence 
of a broader, wavy equatorial wind, as shown for instance in Figure |3] 



Figure |9] suggests that horizontal shear (= barotropic) instabilities are important dy- 
namical ingredients of the sequential flow evolution we have just described. The latitudinal 
profile of zonally- averaged zonal wind velocity, [u], is shown in the top panel, while its second- 
order meridional derivative, d[u]/dy'^ , is shown in the bottom panel, for the same planet day 
5 and a = 0.9 level flow as in Fig. [8j Although necessary and sufficient conditions for 
the development of barotropic (horizontal shear) instabilities are generally not known, the 
Rayleigh-Kuo inflexion point criterion provides a useful necessary condition, which accounts 
for the stabilizing influence of the planetary vorticity ( Kuo|1949 Vallis|2006 ). The instability 



condition is met when d[u]/dy'^ exceeds the planetary parameter {3 = df/dcj) = 2flp cos (j) / Rp 
(latitudinal gradient of the Coriolis parameter), which is shown as a dotted line in the bot- 
tom panel of Figure [9] While the moderate violations of the Rayleigh-Kuo criterion at 
mid-latitudes and beyond may not be very meaningfuQ the strong violations on each side 
of the equatorial wind, at ~ ±20° latitude, are consistent with the substantial zonal shear 
present there and the associated small scale disturbances shown in Fig. [8) 

We note that the significance of barotropic instabilities operating in the flow, as sug- 
gested by Figures. [8] and |9| is that they could play an important dynamical role by tapping 
the free energy available in the horizontal shear flow and thus possibly limit the asymptotic 
speeds of winds in our model. While wind acceleration followed by saturation in the flrst 
~ 10 days, as shown in Figure [5| appears to be broadly consistent with this notion, addi- 
tional flow diagnostics beyond the scope of the present study would be needed to establish 
more confldently this possibility. 



^Away from the equator, the flow is not as strongly zonal as in the equatorial regions shown in Fig. 
Applying a zonal instability criterion in these regions may thus be of limited value. 
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Discussion and Conclusion 



In this work, we have presented a simple Earth-like general circulation model based 
on the IGCM dynamical core. We used this model to contrast the baroclinic circulation 
regime of the Earth's lower atmosphere with the more barotropic circulation regime that 
emerges from a straightforward extension of the model to the atmospheric flow on a hot 
Jupiter. The distinction between these two major (barotropic and baroclinic) regimes of 
atmospheric circulation is an important one in "geophysical" fluid dynamics. For instance, 
both regimes are relevant to the Earth's atmosphere and critical to our understanding of its 
general circulation, with a baroclinic (lower-level) troposphere and a barotropic (higher-level) 
stratosphere. 

Various factors contribute to the degree of baroclinicity/barotropicity of an atmospheric 
flow. The more stably-stratified an atmosphere is (i.e. the more "radiative" it is, as opposed 
to convective, to use the language of stellar physics), the larger its external Rossby deforma- 
tion radius is, the weaker baroclinic instability growth is and thus the more barotropic the 
circulation regime will be (e.g. Pedlosky 1987; Cho et al. 2008). The strong external irradi- 
ation experienced by hot Jupiter atmospheres creates rather strongly-stratified temperature 
profiles in their photospheric regions (e.g. Seager & Sasselov 1998; Sudarsky et al. 2000; Iro 
et al. 2005; Barman et al. 2005; Fortney & Marley 2007). This relative vertical stability, to- 
gether with slow (synchronized) rotation and high atmospheric temperatures, leads to large 



external Rossby deformation radii (Showman & Guillot 2002 Cho et al. 12003 Menou et al. 



2003 Cho et al. 2008) and favors a barotropic circulation regime (with vertically aligned 



horizontal motions in the various atmospheric layers). The lack of surface drag on the at- 
mospheric flow, for otherwise identical forcing conditions, also favors a barotropic regime 



as horizontal shear tends to inhibit the development of baroclinic instabilities (e.g., James 



& Gray 1986 James 1987 Robinson 1997) For all these reasons, a single-layer, vertically- 



integrated barotropic treatment of horizontal motions in hot Jupiter atmospheres may be 



justified (Cho et al. 2003, 2008 Menou et al. 2003 Salby 1989). 



Since the distinction between barotropic and baroclinic regimes depends on the degree 
of atmospheric vertical stratification, which is known for the Earth but a priori unknown for 
remote exoplanets, the results from our shallow hot Jupiter model should only be interpreted 
as suggestive that this regime is relevant to hot Jupiter atmospheric fiows. A more systematic 
exploration of circulation regimes on hot Jupiters will be needed to address this issue more 
thoroughly. The shallow hot Jupiter model presented here is rather specific and idealized 
in a number of important ways. For instance, adopted values for the profile of relaxation 
temperatures and the radiative relaxation time are rather arbitrary. The presence of deeper 
atmospheric layers and their interaction with the modeled layers has also been ignored in 
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this shallow model. 



Nevertheless, the model may capture important dynamical features of the atmospheric 
flow on hot Jupiters. In particular, the simulated flow has a number of similarities with com- 
parable results reported in the literature, together with noticeable differences. As we have 
already emphasized, our simulated flow is characterized by a broad, zonally-perturbed super- 
rotating equatorial wind, large scale polar vortices, unsteadiness and subsonic wind speeds. 
The emergence of a super-rotating (eastward equatorial) wind in this hemispherically-forced 



flow is consistent with the results of Showman & Guillot ( 2002 


); 


Cooper & Showman 


(2005); 


Showman et al. 


(2008 


); 


Dobbs-Dixon & Lin 


(2008 


); 


Langton & Laughlin 


(2008 


). The broad 



width of this equatorial wind and the presence of counter-jets at midlatitudes is also con- 
sistent with the specific results of Showman et al. (2008) and Cho et al. (2003, 2008). As 



discussed in detail by Cho et al. (2008), however, the typically opposite (westward) direction 
of the equatorial wind that emerges in equivalent-barotropic simulations may point to some 
limitation of that approach. On the other hand, the strong zonal disturbances of the equa- 
torial wind in our model is very reminiscent of a similar pattern of large-amplitude Rossby 
waves discussed by Cho et al. (2003; 2008; see also Langton & Laughlin 2007). This is 



qualitatively different from the zonally-symmetric equatorial flow reported by [Showman &: 



Guillot 


( 


2002 


); 



these zonal disturbances of the equatorial wind are closely related to the flow unsteadiness 



observed in our model (see also Langton & Laughlin 2008). 



While our shallow hot Jupiter flow exhibits large-scale polar vortices, their dynamical 



nature is quite distinct from that of the circumpolar vortices discussed by Cho et al. (2003 
2008). The polar vortices shown in Fig. [7| for instance, are anti-cyclonic, not geostrophically- 
balanced and consistently located on the planetary night-side, somewhat eastward of the anti- 
stellar point. This is in contrast with the geostrophically-balanced, cyclonic polar vortices 
that exhibit systematic longitudinal translations in the equivalent-barotropic flows described 



by Cho et al. (2003 2008 ). We interpret this difference as being due to the forcing-dominated 



nature of polar vortices in our hot Jupiter model, which starts at rest, rather than a dy- 



namical origin like in the flows of Cho et al. (2003, 2008), where polar vortices emerge from 



turbulent initial conditions subject to an energy cascade to large scales under the constraint 



of potential vorticity conservation (see Cho et al. 2003 , 2008 for a discussion; see also Lang- 



ton & Laughlin 2007). We note that, even though our simulated flow is clearly unsteady 
(see Fig. [T]), a possible consequence of the different nature of polar vortices in the present 
shallow hot Jupiter model could be a reduced level of disk-integrated variability, from less 
dynamically active polar vortices, by comparison to what equivalent-barotropic results have 



indicated so far ( |Cho et al.||2003t [Menou et al.||2003[ |Cho et al.||2008t [Rauscher et al.||2007 
2008|. 
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An important difference between our results and several comparable studies reported 
in the literature is the consistently subsonic value of wind speeds in our shallow hot Jupiter 



model. In this respect, our results stand out by comparison with those of Showman fc Guillot 



(2002); Cooper & Showman (2005); Showman et al. (2008); Dobbs-Dixon & Lin (2008); 



Langton & Laughlin (2008). It is presently unclear what is the origin of this fundamental 



discrepancy. We simply note here that one element of answer might be the emergence of 
barotropic (horizontal shear) instabilities in our shallow model, which appear to result from 
the acceleration of the equatorial wind and its flank counterjets and could possibly limit the 
asymptotic wind speeds in our model. More work, including model inter-comparisons, is 
needed to clarify this point. 

The magnitude of wind speeds on hot Jupiters is an open problem. Unlike the atmo- 
spheres of terrestrial planets, giant planet atmospheres lack the large-scale sink of energy 
and momentum that is associated with friction on a solid surface]^ As a result, the only 
sources and sinks of energy and momentum in a hot Jupiter flow as simulated here are the 
Newtonian relaxation, which represents large-scale sources and sinks of radiation, and the 
hyperdissipation, which operates on small scales. The absence of a clearly identified large- 
scale sink of energy (ground friction) makes giant planet atmospheres possibly more difficult 
to understand than their terrestrial counterparts. Indeed, dissipation in the flow, which 
ultimately determines asymptotic wind speeds, is then more strongly dependent on the flow 
itself, via its large-scale coupling to radiation and its effective turbulent dissipation on small 



scales. Goodman (2009) has recently suggested that internal friction between atmospheric 



layers could play an important role in the hot Jupiter context. 

The large variety of flow behaviors found so far in distinct hot Jupiter studies suggests 
that, in addition to model validations such as our simple Earth-like model, inter-comparisons 
of hydrodynamic solvers on identical, well-defined atmospheric circulation problems may be 
necessary to build a solid understanding of these new circulation regimes. The shallow hot 
Jupiter model presented here could be used as a first step in this direction. 



^On the Earth and in our simple Earth-hke model, for instance, the atmosphere reaches a global state 
of momentum balance with the bulk planetary rotation through surface drag, with positive and negative 



contributions depending on the easterly or westerly nature of surface winds (e.g.. 



Holton 


1992 




James & 



Gray|[T986l ). 
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Table 1. Model Parameters 



Parameters 


Model 




T^ii^ vf h - 1 1 K p 


l~rof Tnnif pr 


g (gravitational acceleration [m s J) 


9.81 


8 


lip (planetary rotation rate [rad s Jj 


7.292 X 10"^ 


2.1 X 10"^ 


Up (planetary radius [mjj 


6.371 X 10^ 


10^ 


/c (perlcct gas constant [MKbJj 


287 


3779 


K (— TZ/r fMKSl'l 


0.286 


0.286 


Resolution (T-horizontal; L-vertical) 


T42L15 


T42L15 


2Ndel (hyperdissipation order) 


8 


8 


i^diss (hyperdissipation value [m^ s~^]) 


1.18 X 10^^ 


6.28 X 10^^ 


Tfric (Rayleigh friction time - bottom layer [planet days]) 


1 


oo 


Trad (Newtonian relaxation time - all layers [planet days]) 


15 


0.5 


ATep (equator-pole difference for Teq [K]) 


60 


300 


Ttrop (tropospheric lapse rate for Teq [K m~^]) 


6.5 X 10-=^ 


2 X 10-^ 


Tgurf (base value for Teq [K]) 


288 


1600 


-s^stra (height of tropopause for Teq [m]) 


1.2 X 10^ 


2 X 10^ 


STstra. (tropopause temperature increment for Teq [K]) 


2 


10 
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Fig. 1. — Cylindrical maps of temperature and velocity at planet day 150 in the Earth-like 
model with T42L15 resolution. The color scale shows temperatures in K. The top panel 
illustrates the strong baroclinic activity present at mid-latitudes in the bottom model layer, 
at the (J = 0.97 level. The bottom panel illustrates the reduced level of baroclinic activity 
and the formation of two mid-latitudinal jet streams in the upper troposphere, shown here 
at the a = 0.37 level. 
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Fig. 2. — Zonally- averaged contours of zonal wind speed ([U] in m/s, top panel) and temper- 
ature ([T] in K, bottom panel) for the same Earth-like model as shown in Figure [T| at planet 
day 150. Sigma is the pressure normalized to the bottom value (1 bar). These results are in 
qualitative agreement with the Held-Suarez benchmark, even though quantitative differences 
exist (see text for details). 
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Fig. 3. — Cylindrical maps of temperature and velocity at planet day 100 in the shallow 
hot Jupiter model with T42L15 resolution. The color scale shows temperatures in K. Flows 
in the top (cr = 0.7 level) and bottom {a = 0.37 level) panels share strong similarities. 
This vertical flow alignment and the lack of clear baroclinic activity is consistent with a 
barotropic flow regime. Note the presence of large circumpolar vortices. Advection of heat 
away from the central sub-stellar region occurs both westward (at the equator) and eastward 
(at mid-latitudes). 
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Fig. 4. — Zonally-averaged contours of zonal wind speed ([U] in m/s) for the same shallow 
hot Jupiter model as shown in Figure |3j at planet day 100. Sigma is the pressure normalized 
to the bottom value (1 bar). Over the entire vertical extent of the modeled region, the flow 
is characterized by a super-rotating equatorial wind and slower counterjet at mid-latitudes. 
Zonal-average wind speeds are well below corresponding sound speeds (which are ^ 2500 
m/s). 
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Fig. 5. — Time evolution of representative velocities and temperatures in the same shallow 
hot Jupiter model as shown in Figs |3] and |4j All quantities are shown at specific locations 
on the cr = 0.5 level. Flow steady-state is achieved after ~ 20 days. Top panel: zonally- 
averaged (solid line) and maximum (dashed line) value of the zonal velocity u along the 
equator. Bottom panel: temperatures at the sub- and antistellar points (top and bottom 
solid lines, respectively), at the equatorial east and west limbs (top and bottom dotted hues, 
respectively) and at the north and south poles (two dashed lines). Asymptotic wind speeds 
are subsonic. Temperatures are far from being horizontally homogenized. 
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Fig. 6. — A test of numerical convergence. Cylindrical maps of temperature and velocity at 
planet day 100 in shallow hot Jupiter models with T21L5 resolution {a = 0.5 level, top panel) 
and T170L20 resolution {a = 0.52 level, bottom panel). The color scale shows temperatures 
in K. The overall similarity of these two maps indicates good overall convergence. Even the 
T21L5 flow shares many of the global attributes of higher resolution simulated flows (see 
text for details). 





Fig. 7. — The flow unsteadiness is illustrated by comparing two temperature and velocity 
cylindrical maps at planet days 97 and 98 in the shallow hot Jupiter models with T85L20 
resolution (at the a = 0.52 level). The color scale shows temperatures in K. Both the 
temperature and flow fields exhibit significant variability over one planet day. 
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Fig. 8. — Cylindrical map of temperature and velocity at planet day 5 in the shallow hot 
Jupiter model with T42L15 resolution, at the a = 0.9 level. The color scale shows tempera- 
tures in K. The map zooms in a specific equatorial region (±40° in latitude, +20 to +160° in 
longitude), where the accelerating equatorial fiow exhibits strong horizontal shear and small 
scale disturbances. 
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Fig. 9. — Zonally-averaged zonal wind {[u], top panel) and its second derivative with respect 
to latitudinal length /dy^, bottom panel) in the same shallow hot Jupiter model with 
T42L15 resolution as shown in Fig. [s] (planet day 5 and a = 0.9 level). The Rayleigh- 
Kuo necessary condition for barotropic instability is satisfied when d[u\/dy'^ exceeds the 
planetary (3 parameter (shown as a dotted line in the bottom panel). The strong violation at 
±20° latitude is consistent with the substantial horizontal shear and associated small scale 
disturbances seen in Fig. IS} 



